Electron kinetics in isolated mesoscopic rings driven out of equilibrium 
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Kinetic theory of nonlinear current response to an external field is developed for mesoscopic normal 
metal rings threaded by a magnetic flux. General expressions for direct current (DC) are derived for 
a non-equilibrium regime. These expressions describe simultaneously a contribution to DC made by 
a non-equilibrium (external) field, as well as the contribution caused by interaction with intrinsic 
fields. Contributions of electron-electron and electron-phonon interactions to the direct current in a 
non-equilibrium systems are estimated. The kinetic equation for electrons in a pumped disordered 
metal ring is solved with taking into account the Coulomb and electron-phonon interactions. This 
gives an estimate of an overheating of the pumped electron system. 
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I. INTRODUCTION 

Usual objects for experimental and theoretical stud- 
ies of mesoscopic phenomena in linear or non-linear 
conductivity are systems with two or more external 
leads. Because of experimental difhculties, mesoscopic 
systems without leads have received less attention, ex- 
cept for a long-standing puzzle of persistent currents 
(PC) in norjnaLpjpactal rings in the Aharonov-Bohm 
configuratiorou'ttQLllj'D. As the persistent current is an 
equilibrium mesoscopic phenomenon, most of theoretical 
papers deal with the thermodynamical properties of the 
electron system; it is assumed usually that the thermody- 
namic equilibrium has been sustained in the experimental 
investigations cited. 

Although an early theoretical studyB of non- 
equilibrium conductivity in mesoscopic rings has pre- 
dicted a peculiar effect of a direct current (DC), which is 
induced by an external alternating (ac) electromagnetic 
field, there has not been perforrpjed experimental search 
for this effect until very recentlyO. 

The interest to non-equilibrium mesoscopic phenom- 
ena in systems without leacLaJias increased recently in 
connection with a hypothesisO'tj that the observed mag- 
nitudes of PC may include a large DC contribution, in- 
duced by an uiiGontrolled external noise field. It has been 
supposed alsdaO that this uncontrolled noise may be the 
reason of the anomalous electron dephasing rate discov- 
ered in experimental studiesa of weak localization cor- 
rections to the conductance of long wires. Comparing 
the dephasing rate caused by an external electromagnet: 
fielctj in long wires, and the magnitude of DC currenti 
induced in a ring by the external field with the same pa- 
rameters (both quantities are linear in the external field 
intensity), a universal relationshira-between the two quan- 
tities has been proposed recently^. However, the situa- 



tion in this field remains still unclear and requires further 
research, both experimental and theoretical. 

This set of problems together with the arising exper- 
imental interest, necessitate a deeper study of quantum 
kinetical processes in systems without external leads. 
The present paper is devoted to this study. An impor- 
tant physical difference between systems with and with- 
out leads is that in systems with leads, one may consider 
electron transport through the system even neglecting 
processes of inelastic scattering of electrons. The role of 
a "thermostat" may be played by scattering processes of 
electrons within the leads or in the external part of the 
electric circuit. On the contrary, in systems without leads 
(like metal rings) the meaning of a steady-state response 
to an external field is not well defined unless dissipation 
processes are taken explicitly into account. 

In an earlier papeilll we have considered these features 
of small mesoscopic systems in the situation where the 
main channel of the energy dissipation is provided by 
the electron exchange with neighboring reservoirs, so that 
the system allows for a simplified "dynamical" descrip- 
tion in terms of the free electron model. Here we develop 
a self-consistent kinetic theory of nonlinear current re- 
sponse to an external field in mesoscopic normal metal 
rings threaded by a magnetic flux. The theory takes into 
account inelastic electron scattering processes, first of all 
the processes of the electron-phonon interaction, which 
are responsible for establishing a steady-state regime in 
the excited electron system. With the use of the Keldysh 
technique we have derived a generic expression for the 
direct (zero-frequency) current, which includes terms of 
two kinds: "kinetic" terms vanishing in the equilibrium 
and "thermodynamic" terms, which describe the interac- 
tion induced correction to PC. These two groups of terms 
have different analytical structure, that corresponds to 
the physical difference between retarded nature of the 
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response current and equilibrium processes; the latter al- 
low for a weighted combination of emission-absorption 
processes without cajAsality requirement for "absorption 
to precede emission" EJ. 

This approach provides a regularized description of the 
system. Moreover, it allows to resolve problems con- 
nected with anomalous "diffusipa" modes of zero fre- 
quency and momentum (in ref.E3 we call them "loose 
diffusons"): we show that for the distribution function 
obeying the steady-state kinetic equation, the contribu- 
tion of "zero-diffusion" modes vanishes identically (see 
Appendix). However, this modes are important in tran- 
sient regimes. 

In the presence of an external field the difference be- 
tween PC and DC parts of the zero-frequency current 
has a somewhat conventional meaning. This is due to 
the abovementioned strong non-linearity of the electron 
system, which restricts validity of the formal expansion 
into powers of the external field. In fact, the PC part of 
the zero- frequency current, defined formally as propor- 
tional to the zero-th power of the external field, is mod- 
ified as compared to its equilibrium value. This occurs 
due to a field-induced change of the steady-state distribu- 
tion function of electrons. Such an "incoherent" kinetic 
effect is an unavoidable consequence of the external field 
action. 

A steady-state regime in systems without leads may 
only be achieved via a mechanism of energy transfer 
from the pumped electron system to an external thermo- 
stat. As is usually assumed, this transfer is provided by 
electron-phonon interaction. The efficiency of this mech- 
anism decreases with decreasing temperature that may 
lead to an overheating of the electron system (note, that 
more efficient inelastic electron-electron collisions con- 
serve energy of electrons and therefore cannot provide 
cooling of electron system). 

Study of field-induced modifications of the state of 
electron system is the major task of the present pa- 
per. In Sec. II we consider in detail the kinetic equa- 
tion (KE) for electrons excited by an external field and 
cooled via interaction with a phonon bath. We have 
found a solution to this KE in the limit of relatively fre- 
quent inelastic electron-electron collisions. In this adia- 
batic limit the state of electron system may be consid- 
ered as a quasi-equilibrium one, described by the Fermi 
distribution function with an effective temperature. This 
state relaxes slowly to a steady-state regime via electron- 
phonon inelastic scattering processes. We have found an 
explicit expression for the steady state electron temper- 
ature. We analyze also to what extent the absence of a 
global equilibrium between electron and phonon systems 
may modify the value of the direct current. 

The paper has the following structure. General consid- 
eration of equilibrium and non-equilibrium direct current 
is carried out in section II. Model description, notations 
and formulation of the problem are given in the subsec- 
tion II A. Formal expressions for various contributions to 
the direct current are derived by means of the Keldysh 



technique in subsection II B. 

In section III we study a field-induced non-equilibrium 
energy distribution of electrons. Kinetic equation for 
electrons with taking into account electron-electron and 
electron-phonon interaction is described in subsection III 
A. A steady state solution to this equation in the adia- 
batic regime and the electron temperature relaxation rate 
are given in the subsection III B. A field- induced over- 
heating of electrons and its physical consequences are dis- 
cussed in subsection III C. 

In section IV we estimate additional contributions to 
direct current caused by the absence of a complete equi- 
librium between electrons and phonons (subsection IV A) 
and by a deviation of electron distribution function from 
the equilibrium one (subsection IV B). We show that 
the latter contribution vanishes in the constant density 
of states approximation (for systems with electron-hole 
symmetry) while the former one is quite small. These 
study justifies the applicability of a free electron model 
for description of DC induced by an external field of a 
moderate intensity. 

In section V we summarize the obtained results. 



II. EQUILIBRIUM AND NON-EQUILIBRIUM 
DIRECT CURRENT 

A. A model and formulation of the problem 

We consider a normal metal ring of circumference L 
and cross-section area S, threaded by a magnetic flux 
if. The ring is assumed to be thin, i.e. <C L. It is 
assumed also that L ^ I, where I is the electron mean free 
path with respect to elastic scattering by impurities. For 
a time independent magnetic flux ip, the ring is known to 
possess an equilibrium persistent current Ipc- A time- 
dependent magnetic flux would lead to the direct (zero- 
frequency) current along the metal ring, this cuprent is a 
nonlinear response to the external perturbationQ. 

Before specifying a particular mechanism of the exter- 
nal action, we emphasize a basic feature of the response 
problem in systems without external leads. A steady- 
state regime in such systems may be reached only due to 
a balance between an incoming energy flux, pumped to 
the system, and an outgoing energy flux to a surround- 
ing thermostat. Therefore, the mechanism of the energy 
dissipation in systems without leads is vitally important 
and should be explicitly taken into account; the response 
problem becomes a kinetical one. This is in contrast to 
systems with leads where the energy balance may be pro- 
vided by outgoing electrons and the system allows for a 
dynamical description (see discussion inEj). 

Below we shall consider a kinetical problem of the sys- 
tem response to an external field with taking into account 
energy dissipation processes. We assume that the "cool- 
ing" of electrons is provided by interaction with phonons 
which transfer efficiently the excess energy to the sur- 
rounding thermostat (note that the electron-electron in- 
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teraction itself does not change the total energy of the 
electron system) . As to the field that excites the electron 
system, the analysis we present below may be applied for 
various physical mechanisms. 

In particular, in_cpler_tp keep the connection with 
our previous studieg3'E£lE3li3, we shall consider the action 
of an external electromagnetic field with the component 
£{t) = —(l/c){d/dt)A{t) along the circumference. Such 
a field may arise, for instance, if the magnetic flux con- 
tains a time dependent part. We represent the vector 
potential A{t) in the form: 
1*00 

A{t) ^ / exp (-iwt) + c.c] (1) 

"'0 

and assume that this external field has the Gaussian 
statistics with the spectral correlation function Sa{^) de- 
termined by 

< A{u)A*{u') >= SA{i^)S{uj - uj'). (2) 

In case of a narrow spectral width (quasi-monochromatic 
field) the above expressions reduce to A{u!) — Ao6{uj—uja) 
and S'yi(w) = |Aop5(w — ^0)7 respectively. 

We are interested in the response to the perturbation 

V{t) = ~{e/c)A{t) J dr*^(r)v*(r) , (3) 

where \I''l'(r) is the electron annihilation operator and v = 
—iV/m — eA/c is the electron velocity operator (m is the 
electron mass). More specifically, we are interested in a 
field-induced electron current of zero frequency (i.e. a 
direct current) flowing along the ring. 

In addition to this particular mechanism of the exter- 
nal action, we shall analyze also a more general situation, 
assuming only that due to some "external perturbation" 
there arises a non-equilibrium steady-state in electron- 
phonon system of the sample, and we shall study how 
this influences the direct current. A physical origin of the 
perturbation may be also purely "intrinsic" , caused for 
instance by the energy release in course of a slow relax- 
ation of the non-equilibrium nuclear system of the sam- 
ple. In this connection we also mention the slow ortho- 
para conversion of the H2 molecules inside a metal. The 
corresponding heat release of 1 ppm of H2 in Cu at A K 
could be as large as 5 x 10""^ nW/g even weeks after a 
cooling-down. The similar heat release is observed from 
the quartz glasai3. 

B. Direct current, general expressions 

First we shall discuss general properties of the direct 
current in the Aharonov-Bohm geometry. We shall use 
the Keldysh formalism which allows to take simultane- 
ously into account interactions with external and internal 
fields (see, e.g., a reviewllZl). The current flowing along 
the ring is given by 

I{t)^-t^Tr{v,G^{t,t)}, (4) 



where e = — |e| is the electron charge, Vx = —i{yx — 
eA/c)/m is the electron velocity operator (x-axis is cho- 
sen along the circumference), and = (G)i2 is the 
Keldysh component of the Keldysh matrix Green's func- 
tion in the standard triangular representationO, and the 
trace is taken over spin and space coordinates. In a 
steady state of the electron system in the absence of in- 
teraction and time dependent external fields, the Fourier 
transform of G^ {t — t') is given by 

G^iE)^fE[G§iE)-G^iE)] , (5) 

where Gq'"^^ (E) is the retarded (advanced) Green's func- 
tion, and Je = ^ ~ 2nF{E) is connected with the elec- 
tron energy distribution function np{E). If the consid- 
ered steady state corresponds to the thermal equilibrium 
with temperature T, = tanh_E/(2r). An equilibrium 
persistent current for a non-interacting electron system 
is given by 

^° ^ I / IS ^^^^^^ ~ ^° ^ • 

In the presence of an interaction or an external pertur- 
bation, we have instead of Eq. (^) 

G^(i^)= (Go + GoMo + -)^^ (7) 

where E is a self-energy matrix. In the minimal (second) 
order of the perturbation theory, the self-energy splits 
into independent parts £ = S^-l-£e-p/i+Se-ej which cor- 
respond to different physical perturbation mechanisms 
(an external field, electron-phonon and electron-electron 
interactions, respectively). Corresponding contributions 
to the direct current / = I{t) have very similar structure, 
so it is sufficient to write down explicitly only expressions 
for the case of electron-phonon interaction. 

1. Electron-phonon interaction 

The time-independent current may be represented in 
the form: 

I = W)=h+h+h + h- (8) 

Here 

/i = ^/ ^M^. [G^{E)gG\E-u)gG\E)- 

G''{E)gG''iE - io)gG^{E)] {[/e - /£-Ji^^(^) 

+ [fEfE-u.-l]AD{Lu))}, (9) 

where g is an operator vertex of the electron-phonon in- 
teraction ^^'^^I'u that will be specified later, u = {uq} 
is the displacement vector; G^'"* are electron Green's 
functions in the absence of interactions (for simplicity. 
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here and below we omit the corresponding subscript "0" ) ; 
Al) = D^-b^] and D^^ are components of the 

phonon matrix Green's function _D: = ~i < 

[uq,(1), U/3(1')]+ >; etc. In the short-hand notation the 
integrand in Eq.(H) may be represented as Tr{Ji}, where 
Ji has the following structure 

Ji = (RAA - RRA) [(/ - 

+ (//'-l)Ai?], (10) 

The parts I2 and have the same integral representation 
as Eq.(^) but with different trace structures instead of 
Eq.(0): 

J2 = [RARZ?^ -I- ARAL>'^](1 - //') ; (11) 



RRR[/i:)^^' 
AAA[/i:i^ 



(//' - l)i^^] 
(//' - l)D^] . 



(12) 



Below we shall treat the phonon subsystem as a large 
reservoir in equilibrium. In this case 



{lj) = N{uj) D'\lo) - D'^iu) 



(13) 



N{lo) = 1 -|- 2nB{oj), where ns(w) is the boson (phonon) 
occupation number; in the thermal equilibrium the func- 
tion N{uj) = coth[a;/(2T)]. The part h describes 
the contribution to the direct current that, being aver- 
aged over disorder, is not exponentially small at high 
frequencies of bosons ui ^ Ec, where Ec — D / is the 
Thouless energy, D is the diffusion coefficient of electrons. 
This is because the combinations 'R{E — oj)A{E)'R{E) 
and A(i?pj-a;)A(£^)R(£') allow to build a zero-frequency 



cooperonEl. The part I2 (|ll|) is ra, boson analog of the 
Ambegaokar-Eckern contributiontS to the averaged per- 
sistent current that is exponentially small at high fre- 
quencies of bosons. Finally, the part I3 does not con- 
tribute to the disorder averaged current but exhibits only 
in mesoscopic fluctuations. 

A beauty of the expressions is that the "kinetic" part 
Ii is identically zero in the complete equilibrium. Indeed, 
using ( |l3| ) we can represent Eq.(|lO|) in the form: 

Ji = (RAA - RRA) AL* x (14) 
[(/-/')iV(c.) + //'-l], 

where the last square bracket vanishes in the complete 
equilibrium. The existence of the kinetic contribution 
Eq.(|l^ is the main and generic difference between a 
non-equilibrium steady state and the complete thermo- 
dynamic equilibrium. 

2. Electron-electron interaction 

Expressions for /1-/3 are determined by Eqs.(|9|)-([l2|) 
with substitutions: g — > 1 and D.{^) ^ il(^)j where V. is 
the Keldysh matrix for the screened Coulomb interaction, 
= l/o/[l_yo7r«(^)] and V^^ = V^tt^V^- here Vo 
is the usual (unscreened) Coulomb potential and tt is a 
(matrix) polarization operatoiO (see section IV. B). 



3. Interaction with an external field 

The contribution to the direct current made by an ex- 
ternal classical field does not contain a "thermodynami- 
cal" part Eq. (|ri|) . This is due to general (causal) analyt- 
ical properties of the response: there may be only zero 
or one change of the analyticity [R <-> A) between the 
beginning and end points of the electron lines (see e.g.EJ 
and references therein). Expressions for the field- induced 
direct current are given by Eqs.(||), (|lO|), and ( [T^ ) with 
the following modifications: g — > — euj/c, 

b^{uj) ^ -2mSA{^) , (15) 
and omission of terms proportional to (//' — l)!)^^"^^. 

4-. External field vs. non- equilibrium nodes 



One can easily see from Eqs.(|_3|) and (|15|) that the ef- 
fect of an external field is similar to the action of a non- 
equilibrium part of the phonon field described by a non- 
equilibrium component of the phonon energy distribution 
function Nneq{^) = N{uj) — coth [w/(2r)]. We arrive at 
an important conclusion: the absence of a complete equi- 
librium in the electron-phonon system of a mesoscopic 
ring may result in a non-equilibrium direct current simi- 
lar to the DC induced by an external classical field. 

In the equilibrium, the sum of I2 + I3 gives the quan- 
tity CTeg from Ref.M (jifiwe use the relationship D^{uj) — 
-D'^iuj) used in Refllil). 



III. NON-EQUILIBRIUM ENERGY 
DISTRIBUTION IN A PUMPED ELECTRON 
SYSTEM 



A. Kinetic equation 



From Dyson equations ^Gg — G = I and 
G (^Gq — E j = / for the matrix Green's function it fol- 



lows 



Gq , G 



S, G 



(16) 



where Gq (1, 1') = [idt, - H{ri)]6{l - 1')- Taking coor- 
dinate trace from the Keldysh component of the above 
equation we obtain: 



2TTvVdtfE{t) = Tr{AI]G^ - AG} 



(17) 



where v is the averaged electron density of states at the 
Fermi-level, V is the sample volume, AS = YJ^ — TA^ and 
the averaging over disorder realizations is implied on the 
right hand side. Using the second order expression for S 
we arrive at the kinetic equation in the standard form 



dtfsit) = If+Ie-ph + Ie 



(18) 
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Here the first term on the right hand side describes pump- 
ing by the external electromagnetic field 



1 



where 



Jo 



(19) 



(20) 



is a normalized spectral distribution of the external field 
and 

1 



2D{e/cfA^{t)=^D{e/cf S{uj)duj (21) 
T-/ Jo 



is a field induced dephasingH 

The second, electron-phonon collision term in Eq.(p^) 
is given by: 



-■e—ph 



^ / dvdv' 



iNiL0)[fE-fE-^] + fEfE-^-l) (22) 



Here 



where 



[/(k) = (1/V) Uaik) exp i-ikra) , (26) 



Uaik) = J drUair)exp{-ikr) (27) 

is the Fourier transform of the potential of the a-th im- 
purity. For the sake of completeness we write down also 
an expansion of the displacement u(r) over the phonon 
creation and annihilation operators 



(ajiq) exp [i(qr - Wj(g)i)] -|- H.c.) (28) 
and the expression for the phonon Green's function 

AI?„^/3(1,1') = 

^ E / ^) - ^t)] ; (29) 



1. Model of the electron-phonon interaction 

We assume that the mesoscopic ring is embedded into a 
medium with a high thermal conductivity. This allows us 
to consider the phonon subsystem of the ring as an equi- 
librium thermal bath. To simplify the problem, here we 
do not take into account effects of phonon scattering at 
the interface, the presence of surface (interface) phonons, 
etc. The electron interaction with long-wavelength crys- 
tal lattice displacements (acoustic phonons) is treated in 
the "jelly model" and is determined by the standard de- 
formation potential interaction 



lef 

e—ph 



iC 



^c^(p + q)c(p)(quq) 



(23) 



where the deformation potential C — ppVp/3, and c 
are the electron creation and annihilation operators, and 



Uq is the Fourier transform 



Uq = y"drexp(-fqr)u(r). 



(24) 



In systems with electron scattering by impurities (i.e., in 
real metals) one has to take into account also the effoct-pf 
phonon-induced impurity displacements (see Refs.liBEJ). 
This results in an additional part of the electron-phonon 
interaction that does not have free parameters once the 
electron-impurity interaction is chosen 

K'^pH --^Y. Uik)cHp + q + k)c(p)(kuq) . (25) 

P,k,q 



AD«(q,c.) = 

[S{uj - LOjiq)) -d{uj + LO,{q))] T]l^^^{q) 



(30) 



In Eqs.(|28|)-(|3C|) pm is the material density; j accounts for 
the longitudinal (j — I) and two (equivalent) transverse 
phonon polarizations (j — ti and j — t2). For longitu- 
dinal phonons ii^fj{c\} = qaqp/q^] while for transverse 
phonons 



J3) 



Sa.,/3 - qaq/s/q"^ 



(31) 



To perform the disorder averaging of the expression 
Tr{AG{E)gaAG{E-u;)g0} in the collision integral (H), 
one should take into account that the impurity potential 
U enters also the interaction vertex g™^ = —iU{'k)ka, 
that corresponds to Eq. (|5|) . Various contributions to the 
averaged value of Tr{AG{E)gaAG{E — uj)gf3} are repre- 
sented by diagrams shown in Fig.l (they are similar to 
those describing a disorder averaged-electron self-energy 
in the paper by Reizer and SergeevE3). 

First we consider the contribution of transverse acous- 
tic phonons. The deformation part Eq.(23) of the 



electron-phonon interaction is absent in this well 
as the diagrams a, c, e, and f of Fig.l. Moreover, there is 
no effect of "diffuson dressing" of interaction vertices (the 
diagram g): a dressed vertex with an entering phonon 
wave vector q is proportional to q and does not con- 
tribute due to the transverse structure Eq.(^l|) of the 
phonon Green's function. Thus, the leading contribution 
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(a) 



(b) 



(c) 



(d) 




(e) 



(f) 



(g) 



+ 



(h) An effective (renormalized) vertices of the 
electron-plionon-impurity interaction. 




(i) Diffuson proparators. 

FIG. 1: (a-g) - diagrams for the averaged kernel of the 
electron-phonon colhsion integral in the kinetic equations. 
Lines describe electron propagators; dashed and zigzag lines 
correspond to the correlation function of the disorder po- 
tential and to the phonon propagator, respectively. Circles 
and small squares denote vertices of the usual (deformation 
potential) electron-phonon interaction and of the electron- 
phonon-impurity interaction, respectively; large squares de- 
note an effective (renormalized) vertices of the electron- 
phonon-impurity interaction, determined by the diagram h. 



is given by simple ( "Drude" ) diagrams (b and d) and one 
obtains the following expression for the le-ph (22) 



le-ph= / dujK{L0){N{L0)[fE+^+ Je-u. 
<J —oo 

-2fE] + fE[fE+u, - fE-^] + 2) . (32) 



Here the kernel K{uj) is given by 



(33) 



where Vt is the sound velocity of transverse acoustic 
phonons and a dimensionless (transverse) coupling con- 
stant (3t is determined as 

/ Pfvf Y f 

V 3 / p,nV^ 



(34) 



and the function $4 (x) is defined by 

$t(x) 1 + 3[x - {x^ + l)arctan(x)]/(2a;3). (35) 

In the Umiting cases: $t(a;) « at a; ^ 1 and 

$t(a:) « 1 at a: » 1 . The same function arises in stud- 
ies of pJaoHon effects on the conductivity of disordered 
metalSM 

Similarly, one may obtain a contribution of longitudi- 
nal phonons. Now all the diagrams a-g in Fig.l should 
be taken into account. In the limit of small phonon mo- 
menta q -^i l/l, each of the diagrams (e, f, and g) with 
the diffuson installation is considerably larger than the 
corresponding generic diagram (a, c, and d, respectively). 
However, in the leading order in l/(g/)^ the diffuson di- 
agrams mutually cancel and their remaining parts are 
comparable with contributions made by diagrams a-d. 
Such a suppression of the diffuson dressing of electron- 
phononaateraction vertices, demonstrated by Reizer and 
SergeevE3, occurs due to a destructive interference of 
contributions made by the usual deformation potential 
Eq.(23) and by its electron-phonon-impurity counterpart 
Eq.(25): the diagrams e and f are canceled by the cross- 
terms described by the diagram f and its symmetric part- 
ner (not shown) . The physical mechanism of such a can- 
cellation is connected with the Galilean transformation 
from the laboratory system to a system moving together 
with the crystal lattic(£3M 

The resulting expressions for the contribution of longi- 
tudinal phonons differ from Eq.(|33|)-(|35|) by the replace- 
ment (3t —> Pi (where Vt — > vi) and <&t(a;) — + ^i{x), where 



at X <C 1. Typically, Vt is by factor 2-3 smaller 



than vi , that allows one ±0 neglect safely the contribution 
of longitudinal phononsEil. 



B. Steady state solution to kinetic equation 

We are interested in a steady-state regime of the elec- 
tron subsystem subjected to a stationary external pump- 
ing. The only collision channel that takes off the injected 
energy, is the electron-phonon interaction. The electron- 
electron interaction conserves the total energy but leads 
to a redistribution of electrons over energy levels. In a 
closed electron system with a given energy this would ef- 
fectively lead to a quasi-equilibrium Fermi distribution 
with some effective temperature Tg. We shall assume 
that this process of establishing a quasi-thermal distri- 
bution in the electron system is faster than the cooling 
rate due to the electror«honon collisions (similar ap- 
proach was used in RefE3 in studying non-equilibrium 
currents in a metal ring). This allows us to use an adi- 
abatic approach: we put fE{t) ~ tanh[£'/(2Te(i))] and 
N{E) = coth[E/{2T)] into Eq.(|I|) and integrate over 
E with the weight E. Due to the energy conservation, 
the electron-electron collision term makes no contribu- 
tion and we arrive at a closed equation for the effective 
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electron temperature T^. 



(36) 



+ 4 / dLuu^K{uo) 

10 



coth (—] 
\2TJ 



— coth 



2T, 



where ujq is the " mean-square-root" of the external field 
frequency 



,2 _ 



(37) 



In the derivation of Eq.(p6D we have used an identity 
//'-! = -N{E - E'){f - /') for the Bose and Fermi 
distribution functions at the same temperature. 

An effective electron temperature in a steady- 
state regime is determined by the stationary solution to 
Eq. (p6|) . Small deviations ST^ from the stationary solu- 
tion Tg are described by a linearized Eq.(36): 



dSTe/dt = —STe/rx, 
where the temperature relaxation rate is given by 
3 r to^K{Lj) 



1 

TT 



dio . 



(38) 



(39) 



To obtain explicit expressions for Tg and tt-, we restrict 
our consideration to the limiting cases where the electron 
elastic mean-free path I is small or large as compared 
to the phonon wavelength ^ vt/uj which corresponds to 
typical transfer energies to ^ T^. 

Using Eqs.(|3^) and ( |35| ) in the case of small x = uil/vt, 
we arrive at the following expression for the effective 
steady-state electron temperature 



-,1/6 



288C(6)/3tfc|/T/ 



and for the temperature relaxation rate 1/tt'. 

54C(6)/3J(fcsre)4 



1 

TT 



(40) 



(41) 



In these expressions ({x) is the Riemann zeta-function 
(C(6) « 1); for further estimates we have written explic- 
itly the Planck and Boltzmann constants, h and ks- As 
we see from Eq.(40), a noticeable overheating takes place 
at intensities of the external field determined by Eq. ( |2l| ) 
and the condition 



1 3 • WPtKkBTf 



(42) 



At relatively high temperatures T or at a considerable 
overheating (Te ^ T) , the assumption of small ratio x = 
Lol/vt ~ Tgl/vt may be violated. For the opposite case of 



relatively large transfer energy x = tol/vt 
the kernel K(uj) ~ lu and we find: 



T, = 



1/4 



Tel/Vt > 1 



(43) 



and 



TT 



2i6C(4)A(fcBre)- 

TT'^p],Vtl 



(44) 



For the validity of the used adiabatic approach we should 
require the smallness of the temperature relaxation rate 
1/tt as compared to the rate l/re-e of establishing a 
quasi-thermal Fermi distribution due to electron-electron 
collisions. An estimatejJpr Te_e follows from the electron- 
electron collision ternxZl in the considered quasi-lD ge- 
ometry: 



1 



hg 



\J Eel Ec 



(45) 



Here g — Ec/A is the dimensionless conductance (A = 
l/(i^V)) is the mean electron energy level spacing), Ecut 
is a typical low energy cutoff that depends on the rela- 
tionship between Ec and Tg, but at any case EcutT <C 1, 
where t = l/vp is a, mean free time for elastic scattering. 
Comparing the electron "thermalization" rate l/T^-e 
with the temperature relaxation rate 1/tt Eq. ( p| ) (con- 
sidering, for instance, the case of small typical phonon 
momenta ksTe/ (hvt) < i/l), we obtain: 



TT 



hvt 



(46) 



For typical experimental configurations the ring cross sec- 
tion S is comparable with and the ratio T^-el'^T is 
small due to two other factors on the right hand side of 
Eq.(^). This justifies the validity of the adiabatic ap- 
proach. 



C. Numerical estimates 

To get an idea about i^ii overheating effect, consider a 
gold ring. Following ref.Eil, we take the following material 
parameters: (i — 1.4, pp — 1.3 • 10~^^g • cm/s, and Vt = 
1.2 • lO^cm/s. Then Eq.(|40|) may be represented in the 
form 



70mK 



(t^O • 10-9s)l/3 



[(/•105cm-i)(T/ •108s-i)]i/6 



(47) 



where it has been assumed that Tg ^ T. Eq.(p)[) has 
been derived under the assumption of small transferred 
momenta Tejvt < 1/^ that may be represented as 



Te < 90mK/(/ • lO^cm" 



(48) 



If the fraction on the right hand side of Eq . ( p7| ) does not 
exceed unity, the expressions (^^ and ( |4^ ) are compati- 
ble. 
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We see, that the presence of a high frequency external 
pump field may result in a noticeable overheating of the 
electron system. For luq ~ 10^s~^ and / ^ lO^^cni, the 
electron temperature Tg is measured by tens mK even 
for such an external field that leads to a quite moderate 
field-induced dephasing rate on a scale of 10^-10^s~^. We 
note also that Tg depends very weakly (as P^/^) on the 
pumping power P. 

The overheating may put serious obstacles in investi- 
gations of coherent nonlinear phenomena in closed meso- 
scopic samples. 



IV. NON-EQUILIBRIUM EFFECTS ON THE 
DIRECT CURRENT 



As discussed in the preceding section, pumping the 
electron system may result in a noticeable deviation of 
the electron distribution function from the equilibrium 
one. The absence of the complete equilibrium between 
the electron and phonon subsystems may lead to an ad- 
ditional contribution to the averaged direct current. This 
contribution cannot be considered in the framework of a. 
simplified-dynamical treatment used in the earlier worka 
(see alsoE3), it requires the kinetic approach. Here we 
shall study this effect. In the first and second subsections 
we shall consider non-equilibrium phonon and Coulomb 
contributions to disorder averaged DC. Of principal im- 
portance is the non-equilibrium contribution connected 
with the part Ii Eq.(^) of the DC, as this part vanishes in 
the equilibrium. Calculation of disorder averaged quan- 
tities is performed within the usual formalism of weak 
localization theory. The only point to be commented is 
the appearance of a diffuson carrying zero frequency and 
wave vector. This singular diffuson ("loose" diffusonfiJ) 
is a typical feature of non-equilibrium problems in meso- 
scopic systems. Such a diffuson may be built by drawing 
parallel disorder lines embracing the "self-energy" part 
between R and A Green's function (see Fig. 3 in Ap- 
pendix). Thepcple of the anomalous diffuson has been 
studied earlierli3 in the framework of the dynamical con- 
sideration and it is closely connected with the renormal- 
ization of the electron energy distribution function. Here 
(see Appendix) we prove an important identity for the 
anomalous diffuson. Namely, we show that when taking 
into account all the mechanisms contributing to the "self- 
energy" , the diagrams with the anomalous diffuson can- 
cel each other at a steady-state regime. Thus the loose 
diffusons signal on the incorrect or time-dependent en- 
ergy distribution function, and they cancel out as soon as 
the correct steady-state distribution function /e is sub- 
stituted in Eq.(||). 



A. Non-equilibrium electron-phonon interaction 
effects on DC 

The question about a DC induced due to the absence 
of the equilibrium between electrorts-, and phonons, has 
been put forward by Spivak, et a/.c3, who have calcu- 
lated the DC variance. Here we shall consider not the 
variance but the disorder averaged DC (in the presence 
of a magnetic flux pierced the ring). As in the preced- 
ing sections we restrict our consideration to the case of 
transverse acoustic phonons. 

The leading contribution to DC is given by a diagram 
(see Fig. 2) containing a loop of two cooperon propagators 
connecting two "Hikami boxes" (a rhomb and a triangle) . 




FIG. 2: A leading contribution of the electron-phonon or 
electron-electron interaction to a non-equilibrium part of the 
time-independent (DC) current. Wavy lines correspond to 
cooperons (C) or diffusons (D); the zigzag line corresponds to 
the screened electron-electron interaction or the phonon prop- 
agator. In the case of longitudinal phonons both electron- 
phonon and renormalized electron-phonon-impurity vertices 
should be taken into account. For the case of transverse 
phonons the diffuson should be omitted, and the rhomb and 
triangle are directly connected by the phonon propagator (the 
zigzag line) terminated by the renormalized electron-phonon- 
impurity vertices. 

One of the cooperons carries zero frequency and a small 
wave vector k = 27r(n — 2(j))/L along the ring circumfer- 
ence (i.e.,"x"-axis); is a magnetic flux measured in units 
of /ic/|e|. It is the cooperon that provides the dependence 
of the effect on the magnetic ffux. The dependence of the 
second cooperon C(fc — g, to) on the magnetic flux is negli- 
gible at phonon wave vectors q >> 1/L. Expressions for 
the Hikami boxes are given by Fa{q) for the rhomb and 
by kFp{q) for the triangle. The vector function Fa{q) is 
defined (for the RRA succession of Green's functions) as 

Fc.{q) = '5„,,^l^G«(p + q)G^(p) 
p 

+ f E«-^^"G''(?' + 9)[G^(P)]' (49) 
p 
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For both RRA and RAA cases, the corresponding func- 
tions Fa coincide and equal to 



Fa = Oa,x 5 h{ql) 



(50) 



(a longitudinal part ~ qa has been omitted here), where 
, , ^ arctan(x) 3farctan(a;) — xl , , 

H^) = — + y — - (51) 



At large x, h{x) « ■n/{2x). At small x, h{x) « (4/15)x^. 
For m-th harmonic of the phonon induced averaged cur- 
rent we obtain 



(m) _ 



duj Lo^ h? {ijjl / Ut) 



Gtt^ vDut Jo 
[coth(tj/(2T)) - coth(cj/(2re))]ImC(cj/Mt,tj) .(52) 

The range of small momentum transfer q ~ ujl/ut <^ 1 
corresponds to the diffusion regime, where ImC(<;,Ci;) = 
a;/(w^ -I- D^q^). In this case the leading contribution to 
the integral stems from the range u^/{vfI) ^ w <C ut/l 
and we obtain 



15 



(53) 



Therefore we have an estimate (for Te — T > T): 



r(™) 



eE,[{TjE,){TT,){TjEF){T,/ujD)l (54) 



where lujj ^ utpp is the Debye frequency. The factor 
eEc is a typical (mesoscopic) magnitude of a persistent 
current in a single ring, but the expression in the square 
brackets is extremely small. 

Even smaller quantity arises in the case of large mo- 
mentum transfer g = wZ/itt ^ 1 in Eq.(|5^), that may 
take place at relatively large Te, i.e. Te > ut/l, the lat- 
ter quantity is of the order of 100 mK for ut = lO^cm/s 
and I — 10~^cm. In this case we shall use the ballistic 
expression for the "cooperon", that reads: ImC(g, w) = 
LuKvpqY (at bj <^ vpq). As a result, we shall obtain an 
expression which is by a factor [ut/ilT^)]'^ <C 1 smaller 
than Eq.(U). 

We arrive at the conclusion that the non-equilibrium 
contribution of electron-phonon interaction to the flux- 
dependent part of the averaged DC is negligible. This is 
compatible with the smallness of-Lhe photovoltaic current 
variance found by Spivak, et alS3. 



B. Non-equilibrium electron-electron interaction 
effects on DC 



As noted in sec. II, effects of the absence of the complete 
equilibrium in pumped electron systems may result in a 
non-zero contribution Ii to DC. This effect takes place 
for the phonon contribution to Ii studied in the preceding 
section, although the magnitude of the effect turns out to 



be small. Considering the contribution of the electron- 
electron interaction to Ii we deal with the modification 
D,{lo) V.i^) of the (symbolic) expression Ji Eq.(p^) 
(see subsection IIBl): 



(RAA - RRA) [(/ - nv"" 
iff - 1)AV] 



(55) 



Here = Vo/[l - VoTr^^'^'>] and V^'> = V^tt^V^ 

are elements of the Keldysh matrix y_ for the screened 
Coulomb interaction; AV — — V"^; Vq is the usual 
(unscreened) Coulcpib potential; and tt is a (matrix) po- 
larization operatortZI. In dirty metals (in the diffusive 
electron propagation regime) matrix componen t S || Of the 
screened Coulomb "potential" are given by (seeBEl): 



AViq^u;)^-—^, 
Uq'^v 



and 



y^(g,c.) = - 



2iuj 



(56) 



(57) 



where ly is the density of states (at the Fermi energy) 
corresponding to the effective sample dimension, and the 
quantity uJ is defined by 



^ r [I- fEfE-.]dE. 

^ J-oo 



(58) 



To obtain the contribution of the part Eq.(|55|) to the 
disorder averaged DC, we should average the triangle of 
electron Green's functions (G^(E)G^(E - cj)G^(E) - 
G^(E)G^(E-a;)G^(E)) and integrate the resulting 
expression over E. Formally, the leading contribution 
to the kinetic part Ji of the averaged DC would be given 
by the two-cooperon diagram in Fig. 2. However, we will 
show that this contribution, as well as contributions of 
higher order diagrams, vanish after the integration over 
E. Due to the presence of distribution functions with 
the asymptotic properties fp ±1 at i? ^ ±oo, the 
integration over E runs in a close vicinity of the Fermi 
energy. Here we restrict ourselves to the usual approx- 
imation of the weak localization theory, namely, we ne- 
glect a weak dependence of averaged products of elec- 
tron Green's functions on the common part E of their 
arguments; we neglect also the energy dependence of the 
averaged density of states (DOS). In this approximation 
the only £'-dependent part of the integrand is given by 
the square bracket [...] in Eq.(p5|) and we flnd: 



[...]dE ^2LuV^{q,uj) -2ojAV{q,uj) =0. (59) 



(the latter equality follows from Eqs.(p6|) and (pT])). This 
result has been obtained without any assumption about 
the form of the distribution function fp', the derivation 
is based on the identity 



Ue - fE~u]dE = 2u 



(60) 
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valid for an arbitrary distribution function /e (satisfy- 
ing the standard asymptotic conditions). We arrive at 
the conclusion that effects of non-equilibrium and non- 
thermal electron distribution do not contribute to the 
"kinetic" part Ii of the direct current. Of course, there 
is always a "thermodynamic" part I2 of thcj-direct cur- 
rent, considered by Ambegaokar and Eckernla, but this 
part is nonzero even in the complete equilibrium and is 
only weakly changed by a relatively weak pumping. 

It should be emphasized once again that the vanish- 
ing of the Ii part of the electron-electron contribution 
to DC has been obtained in the approximation of the 
energy-independent DOS and the averaged electron tri- 
angle, as well as neglecting disorder-induced correlations 
between the "electron triangle" and the polarization op- 
erator. Results obtained beyond the cpastant DOS ap- 
proximation will be reported elsewhereu. 



complete equilibrium between electron and phonon sub- 
systems. To have a self-consistent description, we have 
examined such non-equilibrium contributions to DC. We 
have found that the phonon contribution to DC caused 
by the difference of phonon (T) and electron (Te) tem- 
peratures, is negligible as compared to the typical PC 
magnitude. We have studied also the contribution to DC 
caused by the Coulomb interaction of non-equilibrium 
electrons, for a no n- Fermi distribution of electrons over 
energy levels. We have found that within the framework 
of usual approximations of the weak localization theory 
(neglecting energy dependence of the electron density of 
states near the Fermi surface), this non-equilibrium dis- 
tribution vanishes identically. This analysis generalizes, 
supports, and shows thie-Jimits of applicability of the re- 
sults of earlier workdallj, where only the field-induced 
contribution to DC was considered. 



V. CONCLUSIONS 

We have considered kinetics of closed mesoscopic sam- 
ples (metal rings) subjected to an external ac pump 
field. To provide a consistent description of excitation 
and relaxation processes in the pumped electron system, 
electron-phonon interaction has to be taken explicitly 
into account. With the use of Keldysh technique we have 
derived general expressions for a time-independent elec- 
tric current along the ring in the presence of a static mag- 
netic flux. These expressions describe simultaneously a 
thermodynamical persistent current (PC) and a direct 
current (DC) induced nonlinearly by the ac field or, in 
general, by an arbitrary perturbation of the equilibrium 
state. The kinetic approach allows to avoid a problem 
of anomalous (loose) diffusongl3. The singularity repre- 
sented by a loose diffuson in the particular case of dif- 
fusive mesoscopic systems, is generic to non-equilibrium 
systems described in the framework of the Keldysh tech- 
nique. The origin of this singularity is the ansatz Eq.(||). 
The singularity arises if the energy distribution function 
Je does not satisfy the kinetic equation (e.g., if it is cho- 
sen to be the equilibrium Fermi distribution). We show 
that the anomalous diagrams with loose diffusons vanish 
in the steady state regime if Je is chosen to make the 
inelastic collision integral vanish. 

Solving the kinetic equation for the electron system, 
pumped by an external field and cooled through the 
phonon bath, we have derived expressions for an effective 
electron temperature Tg. We have found that for appre- 
ciable external field intensities, there may be a noticeable 
overheating of the electron system with Te being on the 
scale of tens mK. This overheating may hinder manifes- 
tations of mesoscopic effects (like PC phenomena). The 
analysis restricts the range of permissible external pump 
intensities. 

In addition to crude thermal effects, there may be more 
subtle phenomena. Namely, the DC may be induced not 
only by the external field but also due to the absence of a 
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VI. APPENDIX: VANISHING OF THE 
ANOMALOUS DIFFUSON CONTRIBUTION 

As mentioned in the text, there m^-be anomalous di- 
agrams with a "loose diffusion" partll3. Such a diffuson 
may be constructed of retarded {G^) and advanced (G^) 
Green's functions in- and outgoing from a self-energy in- 
stallation E, see Fig. 3. 




FIG. 3: A "head" of the anomalous diffuson. 

As this installation conserves energy (in case of an ex- 
ternal field we separate a part with equal numbers of 
absorbed and emitted "photons"), the energy arguments 
of and G'^ coincide, hence the loose diffusion carriers 
a zero frequency. An anomaly arises in the case, where 
the diffusion "head" (G^SG'^) is averaged over the dis- 
order independently of the averaging the rest part of the 
diagram. Due to the restoration of the translational sym- 
metry, the averaged "head" transfers no momentum to 
the loose diffuson. Thus, we arrive at a singularity con- 
nected with a diffuson of zero frequency and wave vector. 
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Our present task is to show that the contribution of 
anomalous diagrams vanishes in a steady state regime. 
Namely, we shall argue for the vanishing of the diffuson 
"head" H{E) 



n{E) Tr [G{E)E{E)GiE)] 



K 

RA ' 



(61) 



where all functions are in Keldysh's space; the trace 
is performed over space coordinates; the subscript RA 
means that we need to keep only the part with outgo- 
ing and incoming G^ lines. Using the relationship 
G^{E) = fEAG{E), we have: 



HiE) = Tr[G'^iE)[-fE^'^iE) 

+ ^'''{E)+fE^^{E)]G^iE)y 



(62) 



Performing here the disorder averaging one should take 
into account only an "irreducible" part (this is marked by 
a superscript irr), i.e., those diagrams which do not in- 
clude an anomalous diffuson ladder built between G^ and 
G^ lines, as this loose diffuson has been already extracted 
from the "head" (see Fig. 3). First, we shall establish a 
useful identity for an "irreducible vertex" JZ^^^ 



TV"{v^,V2;E) ^ Jdr [G«(r, n; £;)G^(r2, r; E)Y 



X63) 



formed by G^ and G^ lines connected with the same site. 
We shall use the following trick. For the moment, we 
treat an infinitesimal positive quantity S in the definition 
of exact electron Green's functions G^''^-' — [E ~ H ± 
iS/2]~^ as a finite constant, with taking the limit S — > -1-0 
afterwards. This allows to deal with a "reducible vertex" 
TZ with no restriction for the arrangement of disorder 
lines on averaged diagrams. This relationship between TZ 
and 7?."''" is shown symbolically in Fig. 4. 





FIG. 4: Relationship between reducible (bold) and irreducible 
vertices. 



The reducible vertex contains an anomalous part with 
a diffusion of zero frequency and momenta. This diffusion 
equals 1/(5. In the leading order in St{E) <C 1, where 
1/t(E) is the rate of elastic scattering, we obtain 



7^^'-''(rl,r2;S) =<5T(S)7^(rl,r2;^). 



(64) 



Note that neither the left, nor the right hand sides of 
Eq.(^) are singular in the limit (5 — > -1-0. Representing 
the Green's functions in terms of exact eigenstates 
and eigenenergies E^ 



G^^^Hr,r';E)^J2 



E~Ef,±i6/2 



(65) 



we obtain the following expression for the reducible ver- 
tex: 

7e(ri,r2;S) = ^ dr G^(r, n; S)G^(r2, r; £;) 

= (V^)AG(r2,ri;i?), (66) 

where AG = G^ - G^. With the use of Eqs.(||) and 
(|66|), we arrive at the final expression for the irreducible 
vertex K^^^ of our interest: 



(67) 



7^''^'■(rl,r2;i?) = zT(i?)AG(r2, ri; i?). 



As neither the left, nor the right hand sides of Eq.(|67 
are singular, we may safely take the limit 5 +0. Using 
Eq.(|^) we find for the diffuson head n{E) 



HiE) = -iT{E)Tr[[fE^''{E) 
- fE^^{E)]AG{E)] . 



j:^{e) 



(68) 



The trace on the right hand side of the last equation co- 
incides identically with that on the right hand side of 
the kinetic equation (^7|). We conclude that the contri- 
bution of anomalous diagrams vanishes in a steady state 
regime, if the correct steady-state distribution function 
Je is used in the ansatz Eq.(||). 
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